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Abstract. We have applied a relativistic Tolman-Bondi 
model of the Virgo cluster to a sample of 183 galaxies with 
measured distances within a radius of 8 degrees from M 87. 
We find that the sample is significantly contaminated by 
background galaxies which lead to too large a cluster mean 
distance if not excluded. The Tolman-Bondi model predic- 
tions, together with the HI deficiency of spiral galaxies, al- 
lows one to identify these background galaxies. One such 
galaxy is clearly identified among the 6 calibrating galax- 
ies with Cepheid distances. As the Tolman-Bondi model 
predicts the expected distance ratio to the Virgo distance, 
this galaxy can still be used to estimate the Virgo distance, 
and the average value over the 6 galaxies is 15.4±0.5 Mpc. 

Well-known background groups of galaxies are clearly 
recovered, together with filaments of galaxies which link 
these groups to the main cluster, and are falling into it. 
No foreground galaxy is clearly detected in our sample. 
Applying the B-band Tully-Fisher method to a sample of 
51 true members of the Virgo cluster according to our 
classification gives a cluster distance of 18.0 ± 1.2 Mpc, 
larger than the mean Cepheid distance. 

Finally, the same model is used to estimate the Virgo 
cluster mass, which is M = 1.2 10 15 Mq within 8 degrees 
from the cluster center (2.2 Mpc radius), and amounts to 
1.7 virial mass. 

Key words: Galaxies: clusters: individual: Virgo - Galax- 
ies: distances and redshifts - Galaxies: elliptical and lentic- 
ular - Galaxies: fundamental parameters - Galaxies: spiral 
- Cosmology: distance scale 



1. Introduction 



In 1990, Fouque et al. (199C) derived an unbiased dis- 
tance to the Virgo cluster, based on a complete sample 
of 178 spiral galaxies, using the B-band Tully-Fisher re- 
lation (Tully & Fisher |1977p . A previous similar study 
based on 110 spiral galaxies (Sa - Sm) was published by 
Kraan-Korteweg et al. ( |1988| ). Soon after, Teerikorpi et al. 



(1992, hereafter T92) suggested that this distance deter- 



mination may have been contaminated by the inclusion of 
background galaxies into the complete sample, although 
we a priori excluded galaxies generally attributed to the 
background M group and W cloud. 

The study of the s truct ure of the Virgo cluster starts 
with de Vaucouleurs ( 1961 ), who identifies the southern 
extension (Virgo cloud X), and the wing (Virgo cloud 
W). This study is extended in de Vaucouleurs & de Vau- 



couleurs (|1973j) who separate the Virgo I cluster (Virgo 
E, S and S') from the Virgo II cloud complex (com- 
posed of Virgo V, X and Y) and the background W cloud 
(comp osed of Virgo Wa, Wb and W sub-groups). Pa- 
turel (1979) applies taxonomy to disentangle these various 
components. Then, Tully ( 1982j ) discusses the separation 
between the Virgo cluster and its southern extension in 
the frame of the Local Supercluster. Later, Ftaclas et al. 



(1984) identify the M group and the N group. Pierce & 
Tully ( h988 [) use the Tully & Shaya infall model (Tully & 
Shaya 1984) to study the velocity - distance diagram of 
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bright Virgo spiral galaxies. Binggeli et al. (1993) clarify 
the situation by listing galaxies from the VCC catalogue 
( 1985[ hereafter VCC) belonging to each of these clouds 
(W, W, M and southern extension). T92 list several re- 
gions in the veloc ity - d istance diagram: A, B, CI, C2, and 
D. Yasuda et al. (1997) also study the 3-D structure of the 
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cluster. Finally, Gavazzi et al. (1999, hereafter G99) use 
distance determinations to identify new groups, named B, 
E, N, and S. 

Please note that a clarification of the nomenclature is 
highly desirable, as the same letters refer to totally dif- 
ferent groups: A can mean the main Virgo cluster or re- 
fer to galaxies with high velocities in front of Virgo and 
falling into the cluster (T92); B can mean the concentra- 
tion around M 49, or a background group in the same 
region (G99), or even a foreground expanding component 
(T92); E can refer to the elliptical component of the Virgo 
cluster or an eastern group (G99); N can be a group iden- 
tified by Ftaclas et al. (1984) or a northern group in G99; 
S refers to the spiral component of the Virgo cluster or a 
southern group in G99; X can mean that the galaxy lies 
within the X-ray contours (Federspiel et al. [1998 , here- 
after F98) or belongs to the Virgo X cloud (de Vaucouleurs 
196l]) ! 



To investigate the structure of the Virgo cluster and 
determine the mean distance of the main cluster, we have 
used the Tolman-Bondi model of the cluster defined in 



Ekholm et al. (1999, hereafter E99) in the spirit of T92, 



for spiral galaxies whose distances were given b y the Tully- 
Fisher relation in _B-band by Ekholm et al. (200C, here- 
after E00) or in iJ-band by G99, but adding the infor- 
mation given by the HI deficiency of Virgo spiral galax- 
ies. We also used the same model for early-type galax- 
ies whose distances were determined by the fundamental 
plane method in G99, or the recent Tonry's compilation 
of Su rface Brightness Fluctuations distances (Tonry et al. 
2001| , hereafter T01). 

Section 2 defines the observable parameters we will 
use for this study. Section 3 introduces the Tolman-Bondi 
model and lists its adopted parameters. Section 4 investi- 
gates how the discrepant Cepheid distance to NGC 4639 
is well explained by our model, and computes the average 
Cepheid distance to Virgo. Section 5 applies our model 
to the determination of the Virgo cluster mass, and com- 
pares it to other determinations. Section 6 investigates 
the structure of the Virgo cluster and lists galaxies which 
do not belong to the main cluster according to our model. 
Section 7 compares the average distances to Virgo adopted 
in the various references we used to what our model pre- 
dicts. Finally, Sect. 8 summarizes the main results of this 
paper. 

2. Observable parameters 

The selected sample has been extracted from the LEDA 
database. It consists of 584 galaxies within a radius of 8 de- 
grees (2.25 Mpc at 16 Mpc) from M 87 (a 20 oo = 12.51381, 
^2000 = 12.3900) and known recession velocities smaller 
than 3000 km/s. Unfortunately, not all these galaxies have 
distance estimates, and we therefore reduced the working- 
sample to 125 late-type and 67 early-type galaxies, with 
distance measurements from G99, E00, F98 and T01. Let 



Fig. 1. Histogram of the recession velocities, in the Lo- 
cal Group reference frame, for the 584 galaxies within 8 
degrees from M 87. The dotted histogram corresponds to 
galaxies with distance measurements 



us describe in more detail some of the observable param- 
eters collected for these galaxies. 

2.1. Coordinates and recession velocity 

These parameters are directly extracted from the LEDA 
database. Coordinates are given in the J2000.0 equinox, 
and helocientric velocities are corrected to the Local 
Group centroid using the Yahil et al. ( [1977 ) formula. A 
histogram of the corrected velocities for the full sample 
of 584 galaxies is presented in Fig. [l]. It clearly shows 
more galaxies on the high- velocity side than expected from 
a gaussian distribution. A similar but smaller effect also 
seems apparent on the low- velocity side. For galaxies with 
distance measurements, the same features appear. 



2.2. Morphological types 

An accurate estimate of morphological types is important 
at least for two reasons: the HI deficiency parameter (see 
below) is calculated by comparing the HI content of a 
galaxy to the average HI content of an isolated galaxy of 
the same morphological type. Any error in the type there- 
fore translates into a corresponding error in the HI defi- 
ciency estimate. Similarly, distances determined by E00 
use the calibration of Theureau et al. ( 1997 ), where the 
intercept of the Tully-Fishcr relation varies with the type. 
Any error in the type changes the galaxy distance value 
accordingly. 

We used three sources of m orpho logical type determi- 
nations: de Vaucouleurs et al. ( 1991 , her eafter RC3), the 
VCC catalogue and van den Bergh et al. ( 1990 ). The qual- 
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ity of RC3 types depends upon its source, coded according 
to Table 3 of Vol. 1 of the catalogue. If the morphological 
type has been measured on large reflector plates (code R) , 
it is adopted; if it comes from PSS prints (codes P or U), 
we adopt an eye average with the two other sources. The 
adopted morphological types differ therefore from those 
given in LED A and used by E00. 

We adopted a numeric morphological type according to 
the RC3 coding scheme. For distance measurements, we 
decided to restrict application of the fundamental plane 
and SBF methods to types between -5 and -1 (E to L+) 
and of the Tully-Fisher method to types between 1 and 
9 (Sa to Sm). Therefore, 7 galaxies were excluded from 
the late-type sample (types -1, and 10), and 2 galaxies 
were excluded from the early-type sample (NGC 4440 and 
NGC 4531, both of type 1). The final sample thus contains 
118 spiral, 43 lenticular and 22 elliptical galaxies. 

2.3. Distances 

2.3.1. Early- type galaxies 

For early- type galaxies, we have used the G99 compila- 
tion of distances, obtained from the fundamental plane 
method in if-band. According to these authors, their ac- 
curacy should be about 21%. These distances are based on 
an assumed Virgo distance of 16 Mpc. To this compilation 
of 55 elliptical and lenticular galaxies, we added SBF dis- 
tances from T01. It contains 35 early- type galaxies in the 
Virgo cluster area, of which 25 are common with G99. 
Tonry's distances are calibrated independently of any as- 
sumed Virgo distance, and the mean shift compared to 
G99 amounts to: 



(d T -d Gg9 ) =0.64 ±0.84 



(1) 



after rejection of NGC 4638 (d G99 = 8.87 Mpc and d T = 
21.68 Mpc). As the mean difference of the two systems is 
not significant, we adopt an unweighted average of the two 
measurements as our distance measurement for this set of 
25 early-type galaxies, and the uncorrected only available 
distance for the 40 other galaxies (30 from G99, 10 from 
Tonry). Please note that both references include more un- 
certain measurements: in G99, they are identified as be- 
ing due to velocity dispersions smaller than 100 km s _1 ; in 
T01, they appear in a separate table, named "table. poor" . 
We have kept these galaxies in our sample, but we gener- 
ally exclude them from our statistical calculations. 

2.3.2. Late- type galaxies 

For late-type galaxies, we also started from the G99 com- 
pilation (59 galaxies retained), obtained from the H-band 
Tully-Fisher method, with a claimed accuracy of about 
16%. We complemented this list with galaxies with B- 
band Tully-Fisher distances from E00 (41) and F98 (109). 



We do not use galaxies from these lists outside of 8 de- 
grees from M 87 or with morphological types outside of 
our adopted range (1 - 9). 

To convert E00 distances into the G99 system, two 
corrections are done: first, we correct the distances to our 
adopted morphological type, in the 20 cases where it dif- 
fers from the LEDA type adopted in E00. Indeed, the 
intercept b of the B-band Tully-Fisher relation derived by 



Thcureau et al. (1997) and adopted in E00 depends upon 
the morphological type. The correction is given by: 



d(T) = d Em X l 0-2x(b(T LEDA )- 6 (T)) 



(2) 



Although the morphological type differences never ex- 
ceed one unit, the correction can amount to 25% in dis- 
tance. 

Addit ionall y, E00 distances are based on the Theureau 
et al. ( 1997 ) calibration, which leads to H Q « 
55 km s -1 Mpc; as E00 adopt a cosmic velocity of Virgo 
of 1200 km s _1 , this implies a Virgo distance of 21.8 Mpc. 
Therefore, the conversion factor to the G99 assumed Virgo 
distance of 16 Mpc is: 



icor 

"eoo 



0.733 d E oo 



(3) 



To convert F98 distances into the G99 system, we sim- 
ply apply the calibration correction derived from different 
adopted Virgo distances, namely 21.5 Mpc and 16 Mpc: 



d c ¥ ° 9 \ = 0.744 d F98 



(4) 



system. 


We find: 




(^EOO — 




1.21 ±0.29 


(^F98 ~~ 


^mean} 


-1.89 ±0.18 


{dG99 — 




1.13 ±0.25 



As we have 65 galaxies with more than one distance es- 
timate, we can make an intercomparison among the three 
sources of distance data, after their reduction to the mean 



(5) 
(6) 
(7) 



where ci mGan is the mean value of available distance 
measurements for a given galaxy. 

All these shifts are significant. This means that our 
first guess of the conversion factors to the G99 system 
was not as successful as we could hope. We therefore re- 
peated it with new conversion factors, until we reach on 
average negligible shifts of the corrected systems to the 
mean one. This leads to the following adopted conversion 
factors, after rejection of two galaxies with discrepant dis- 
tance measurements (NGC 4180 and NGC 4591): 



icor 
"EOO 

.cor 
U F98 
jcor 
"G99 



0.678 d E00 
0.832 d F98 
0.930 d G99 



(8) 
(9) 
(10) 
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We adopt these new conversion factors to put all three 
references into the same system. Then, considering that 
the advantages of the f/-band (low extinction correction) 
are compensated by the smaller slope of the Tully-Fisher 
relation in B-band, we do not weigh each measurement 
and adopt a straight average of the available distance es- 
timates for each galaxy. We will only use these mean dis- 
tances to classify the galaxies among the various classes 
described in Sect. 6. We display in Fig. || a histogram of 
these 118 average distances. The mean value of 20 Mpc 
and the skewness of the histogram clearly show the back- 
ground contamination of our sample. By comparison, the 
histogram of distances for early-type galaxies appears to 
be less contaminated by background galaxies (smaller 
mean distance and smaller skewness). 

2.4- HI deficiency 

As we are concerned with the distance determination of 
individual galaxies in the Virgo cluster, we quantify HI 
deficiency by means of a distance-independent parameter 
DEF, based on the difference between the expected and 
observed logarithm of the mean (hybrid) HI surface den- 
sity, Shi, that is 



DEF=<log£ HI (T obs ))-log£' 



HI 



(11) 



with Shi = ^Hi/c&opt> and where -Phi represents the cor- 
rected HI flux density integrated over the profile width in 
units of Jykms -1 and a op t represents the apparent opti- 
cal diameter in arcmin. For non-detections, an upper limit 
to HI content is estimated by assuming that the emis- 
sion profile is rectangular with an amplitude 1.5 times the 
rms noise and width equal to that expected for a galaxy 
of the same morphological type and luminosity, properly 
corrected for r cdshif t broadening and viewing incli natio n 
(Solancs et al. 1996 ; see also Haynes & Giovanelli 1984 ). 
The adopted standard values of (log Shi) per type are: 
0.24 for Sa, Sab; 0.38 for Sb; 0.40 for Sbc; 0.34 for Sc; and 
0.42 for later spiral types. 

Among our sample of 118 spiral galaxies, we have HI 
deficiency estimates for 106 galaxies, or 90% of the sample. 

It has been suggested that Tully-Fisher distances may 
be underestimated for highly HI deficient galaxies (T92; 
Fukugita et al. 1993). For this reason, we have investi- 
gated possible differences in the rotation velocities of HI 
deficient and HI normal galaxies in the 12 HI deficient 
clusters identified by Solanes et al. ( |200l| ), as well as the 
influence of the HI content on the Tully-Fisher relation- 
ship of the galaxies in our sample. The results obtained 
allow us to conclude that HI deficiency does not affect our 
Tully-Fisher distances. 



3. Tolman-Bondi model 

We use the Tolman-Bondi model of the Virgo cluster 
as defined in E99. Let us recall that the Tolman-Bondi 



model gives an analytical solution to Einstein's field 
equations for a spherically symmetric pressure-free den- 
sity excess, embedded in an otherwise homogeneous uni- 
verse. The parameters of the model are: the observed 
Virgo cluster velocity in the Local Group reference frame 
(980 kms -1 ), the Virgocentric infall velocity of the Local 
Group (220 kms" 1 ), the Virgocentric density profile slope 
(a = 2.85), and the deceleration parameter of the back- 
ground homogeneous universe, taken as an Einstein - de 
Sitter universe for simplicity (q = 0.5). Although there 
is a general agreement about the values of 3 of these pa- 
rameters, the Virgocentric density profile slope is not well 
known, as it should reflect the distribution of the mass 
around the cluster center, not only the galaxy (light) dis- 
tribution. In E99, it has been constrained using 32 galaxies 
whose distances were known using their Cepheids, gener- 
ally measured with the HST, through the PL-relation. For 
details about the model, the adjustment of its parameters 
and the influence of changing these values, the reader is 
referred to T92, E99 and E00. 

For each galaxy, the angular distance to the Virgo cen- 
ter of mass (taken to be the position of M 87) completely 
defines the exact shape of the Tolman-Bondi S-curve in 
the velocity distance diagram, for a given choice of model 
parameters. Then, the observed recession velocity of the 
galaxy (in the Local Group reference frame) gives one to 
three possible distance ratios of the galaxy to Virgo (see 
Fig. U for examples). The first and the third values cor- 
respond to a galaxy falling into the cluster from in front 
or behind. The second value corresponds to a true clus- 
ter member. Assuming that the Tolman-Bondi model can 
explain most of the observed positions of the galaxies in 
the velocity-distance diagram means that we neglect other 
possible components, such as an expanding one suggested 
by T92, or galaxies projected by chance onto the clus- 
ter, but taking part in the Hubble flow. Certainly, not all 
galaxies in the Virgo cluster should follow the Tolman- 
Bondi model. For instance, galaxies belonging to the viri- 
alized core of the cluster do not exhibit a velocity - dis- 
tance relation. It is however difficult to disentangle those 
galaxies from our class 2, so that we include them into it. 
On the other hand, not all galaxies in the class 2 (defined 
as the descending branch of the S-shape Tolman-Bondi 
curve) belong to Virgo: for instance, at large angular dis- 
tances from M 87 the descending branch still exists, but is 
no longer related to proper Virgo galaxies. We assume that 
within 8 degrees of the cluster center, most of the galax- 
ies on the descending branch are true cluster members. 
We will test this strong assumption in Sect. 6 by build- 
ing a histogram of "Virgo distances" , where each point is 
computed from a galaxy distance and its distance ratio to 
Virgo, as given by the Tolman-Bondi model. 

To do this, we must select the most plausible one 
among the possible distance ratios for a given galaxy, when 
the model gives more than one value. Here, we made use 
of two criteria: we try to choose the value closest to the 
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Fig. 2. Histogram of the adopted distances for the 118 spiral galaxies (left) and the 65 early-type galaxies (right) 



assumed Virgo distance (16 Mpc), and, for spiral galax- 
ies, we take into account the HI deficiency, assuming that 
a true member of the cluster has a higher probability to 
be deficient than a galaxy falling for the first time into 
the cluster. We therefore attribute to each galaxy a "class 
number" , which is 1 if we adopt the first value of the dis- 
tance ratio to Virgo (galaxies falling from in front into the 
cluster), 2 if we adopt the second value (true cluster mem- 
bers), and 3 if we adopt the third value (galaxies falling 
from behind into the cluster). 



4. Distance to the Virgo cluster from Cepheids 

There are 6 galaxies in or close to the Virgo cluster whose 
distance is known thanks to the HST observations in V 
and / bands, using the Cepheid period- luminosity re- 
lation. We use the distances published in Freedman et 



al. ( 2001 ) (without the uncertain metallicity correction). 
Five of the six distances cluster about a mean distance 
of 14.6 Mpc, while the sixth one (NGC 4639) gives a 
larger distance (21 Mpc). This is perfectly explained by 
the Tolman-Bondi model, and this was in fact the initial 
motivation of this study. Figure || gives the result of the 
application of the Tolman-Bondi model to these galaxies, 
with a dotted line showing the recession velocity of the 
galaxy, and an arrow showing the distance ratio to Virgo, 
for an adopted distance to Virgo of 15.4 Mpc (see below). 

Two galaxies belong to the Virgo southern extension 
(NGC 4496A and NGC 4536) and lie at more than 8 de- 
grees from M 87. However, their position in the velocity- 
distance diagram should be explained by our model, if we 
assume that no other cluster perturbes them. However, we 
find that the maximum velocity explained by our model 
for NGC 4536 is 1563 kms -1 at 0.83 (iyirgo, too small 



compared to its observed velocity of 1641 kms -1 . For this 
galaxy, our model predicts it to be in the class 3, with 
a distance of 1.79 dyi Tgo . This in turn would lead to an 
inacceptable Virgo distance of dvirgo = 8.07 Mpc, using 
the observed galaxy distance of 14.45 Mpc. In fact, it is 
well known that random velocities about 80 kms -1 exist 
for all galaxies and explain the velocity dispersion of small 
groups of galaxies (Gourgoulhon et al. 1992). We therefore 
assume that this is the origin of the small discrepancy ob- 
served for NGC 4536. 

Table [| gives the adopted distances to the six galax- 
ies, the model predicted distance ratios to the Virgo dis- 
tance for the three possible solutions, the corresponding 
predicted "Virgo distances", the adopted class, and the 
adopted Virgo distance for each galaxy. Now we discuss 
each galaxy in some detail: for NGC 4548, only the second 
and third distances are plausible; but the large HI defi- 
ciency leads us to attribute this galaxy to class 2 as a true 
member of the Virgo cluster. Similarly, the large HI defi- 
ciency of NGC 4321 makes us attribute it to class 2. On the 
contrary, the insignificant deficiency of NGC 4639 leads to 
put it into class 3, with a derived Virgo distance interme- 
diate between NGC 4321 and NGC 4548. For NGC 4535, 
the HI deficiency is low, but putting it into class 1 would 
lead to a larger derived Virgo distance than the previous 
galaxies, so we attribute it to class 2. NGC 4496A belongs 
to the Virgo southern extension and is not HI deficient; 
it is therefore tempting to put it into class 1, but again 
putting it into class 2 leads to a derived Virgo distance 
in better agreement with the remaining galaxies. Finally, 
NGC 4536 also belongs to the Virgo southern extension 
and has a small HI deficiency, and it may be put into class 
1 or 2, once its recession velocity is corrected for a random 
component. The resulting mean value of the six derived 
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Virgo distances is 15.4 ± 0.5 Mpc, which we adopt as our 
estimate of the cluster distance. 

At the referee's request, we have investigated what 
happens if a Virgo distance of 21.5 Mpc is preferred from 
independent arguments. Then, four of the five galaxies 
classified in class 2 now fall into class 1 (infalling galax- 
ies), NGC 4639 is the only true member of the Virgo clus- 
ter among the six, and NGC 4548 position cannot be ex- 
plained by the model. The high HI deficiency of NGC 4321 
is contradictory. This alternative distance is clearly less 
probable according to these six galaxies with Cepheid dis- 
tances. A definitive answer will only become available with 
a larger sample of Virgo galaxies with accurate distances. 



about (1 — 3) x 10 Mq. These values are confirmed by 



Schindler et al. ( |1999| ) who derive 2.1 x 10 14 M within 
1.5 Mpc around M 87 (at 20 Mpc) and 0.87 x 10 14 M 
within 0.75 Mpc around M 49. Our Tolman-Bondi mass 
for the same distance would be M = 1.5 10 15 Mq, almost 
one order of magnitude larger. We do not have any ex- 
planation to offer for this discrepancy, but we note that 
our large mass estimate and steep density profile are sup- 



ported by Tully & Shaya (1998), who find a Virgo mass 
of M — 1.3 10 15 Mq from a modeling of the velocity field 
of the Local Supercluster, assuming a mass-to-light ratio 
of 15OM /L in the field, but 1000 for the Virgo cluster, 
and n o = 0.3. [] 



5. Mass of the Virgo cluster 

Derivation of the Virgo cluster mass from the Tolman- 
Bondi model, and its comparison with the virial mass es- 
timate, have been discussed in T92, E99 and E00. Here, 
we give a summary of the useful formulae and derive the 
Virgo mass for our adopted model parameters. Let i?viR 
be the Virgo distance in Mpc, Vyf™ the cosmic recession 
velocity of Virgo in kms -1 and d the radius normalized 
to the Virgo distance. The mass enclosed within d is the 
product of the "Einstein - de Sitter mass" within the same 
radius by the mass excess due to the cluster. It is given in 
solar mass units by: 



M(d) = M(d) EdS x {l + k'd- a ) 
M(d) EdS = 2.325 10 s q a R ym (V^) 2 d 3 



(12) 
(13) 



k' , the mass excess within d = 1, i.e. at the Local 
Group location, only depends on the observed Virgo clus- 
ter velocity, the Virgocentric infall velocity of the Local 
Group and the cosmological parameters of the adopted 
background universe (q and H Q xt ). It does not depend 
on the Virgocentric density profile slope a. Its value is the 
same as in E99, namely k' — 0.606. For Virgo, we define 
the radius d as corresponding to 8 degrees, which gives 
d = 0.141, in place of 0.105, corresponding to 6 degrees in 
E99. 

With V^f" = 1200 kms- 1 and i? V m = 15.4 Mpc, we 
get M = 1.2 10 15 M^. T he virial mass of Virgo as given 



by Tully & Shaya (1984) is 



M virial = 2.325 10 8 x {nRn) x (3 4) 



(14) 



where notations are those from Tully & Shaya. Their result 
transposed to 15.4 Mpc gives M v i ria i = 6.9 10 14 Mq, so 
that M = 1.7 Mviriai. E99 found a coefficient 1.62 for their 
adopted parameters of the Tolman-Bondi model. 



By comparison, Bohringer et al. (1994) estimate the 



mass of the M 87 sub-cluster from X-ray emission mea- 
sured by ROSAT to - (1.5 - 6) x 10 14 M within a 
radius of 1.8 Mpc at 20 Mpc (5 degrees). Contributions 
from the M 49 and M 86 sub-clusters are negligible, at 



6. Structure of the Virgo cluster 

Disentangling the different components of the Virgo clus- 
ter is difficult and somewhat subjective. However, the 
use of the Tolman-Bondi model allows us to classify each 
galaxy into one of the three classes defined previously. We 
did that independently for early-type and spiral galaxies. 
In each case, we have built several diagrams: the first one 
is the histogram of the "Virgo derived distances" , and it is 
shown in Fig. ||. Both histograms now have a well-defined 
peak around the assumed Virgo distance and are symmet- 
ric, a posteriori confirming our hypothesis. Mean distances 
in both cases are compatible with the assumed distance 
of 16 Mpc. However, the rms dispersion is slightly lower 
in the case of early-type galaxies: this may reveal that the 
spiral sample class attribution is still not fully satisfying, 
or that our restrictive hypothesis rejecting expanding com- 
ponents or projected galaxies in the Tolman-Bondi model 
are not completely fulfilled. 

The second diagram gives the velocity repartition 
among the different classes, and is displayed in Fig. |[ 
In fact, we have split class 3 into two classes, because it 
clearly contains high velocity objects (class 4, above 1600 
km s -1 ), which are background galaxies, and low veloc- 
ity objects (class 3, - 1300 km s _1 ), which are infalling 
galaxies from behind the cluster. We have not convinc- 
ingly identified any spiral galaxy belonging to class 1; the 
only two early- type objects classified in class 1 exhibit low 
velocities (about 800 km s _1 ), which is difficult to inter- 
pret; in fact, these two galaxies may be genuine members 
of the cluster (class 2) with bad distance measurements; 
indeed, both galaxies have a velocity dispersion smaller 
than 100 km s , which makes their fundamental plane 
distances quite uncertain. 

A third diagram gives the repartition of HI deficiency 
among the classes, obviously only for spiral galaxies. It 
is displayed in Fig. ^. It is clear that there is no highly 
deficient galaxies among classes 1, 3 and 4, as expected 



1 Although Tully & Shaya do not specify their adopted dis- 
tance to Virgo, we assume it is the same as in Tully & Shaya 
(|l984|), namely 16.8 Mpc. 
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Fig. 3. Tolman-Bondi model for the six galaxies with Cepheid distances. The dotted line corresponds to the recession 
velocity of the galaxy, while the arrows mark the expected distance ratio if i?vm = 15.4 Mpc (solid arrow) or 
-Rvm = 21.5 Mpc (dotted arrow) 

Table 1. Adopted Cepheid distances to the 6 Virgo galaxies, ratio to the 
Virgo distance given by the Tolman-Bondi model and corresponding Virgo 
distances for each possible class, adopted class and Virgo distance, and HI 
deficiency 



Name 


d 
Mpc 


di 


d/dy 
d 2 


d 3 


class 


dy 

Mpc 


def 


NGC 4321 


14.32 


0.68 


0.95 


1.71 


2 


15.07 


0.49 






21.06 


15.07 


8.37 








NGC 4496 A 


14.52 


0.75 


0.92 


1.76 


2(1) 


15.78 


-0.09 






19.36 


15.78 


8.25 








NGC 4535 


14.79 


0.77 


0.96 


1.91 


2(1) 


15.41 


0.19 






19.21 


15.41 


7.74 








NGC 4536 


14.45 


(0.83) 


(0.83) 


1.79 


1 or 2 


(17.41) 


0.25 






(17.41) 


(17.41) 


8.07 








NGC 4548 


15.00 


0.23 


1.07 


1.30 


2 


14.02 


0.83 






65.22 


14.02 


11.54 








NGC 4639 


20.99 


0.47 


1.03 


1.45 


3 


14.48 


0.10 






44.66 


20.38 


14.48 
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Table 3. List of galaxies out- 
side the main Virgo cluster and 
not attributed to any group, with 
their distance ratio to Virgo 



Name 


d/d v 


Class 


NGC 4417 


0.39 


1 


NGC 4488 


0.46 


1 


NGC 4591 


2.21 


4 


NGC 4598 


1.92 


4 


NGC 4698 


1.44 


3 


IC 3298 


2.27 


4 


IC 3483 


1.22 


3 


UGC 7697 


2.34 


4 



Fig. 6. Distribution of HI deficiency among the different 
classes for spiral galaxies 

because we used HI deficiency as one of our classification 
criteria. 

Finally, Fig. |7] displays the repartition on the sky of 
spirals and early-type galaxies with different symbols for 
different classes. If we compare the distribution of the 
two types of galaxies, it clearly appears that the different 
classes concentrate into common regions. One can easily 
recognize the W group around (12. 3h, 6°) and the M group 
around (12. 2h, 13°). Both groups have a mean distance 
ratio to Virgo of 2.1. Galaxies at the same distance but 
slightly displaced from the main concentration have been 
classified as group halo. Finally, galaxies with intermedi- 
ate distances between those of these groups and the Virgo 
distance, and classified into our class 3, are interpreted as 
belonging to extended filaments extracted from the back- 
ground group and falling into Virgo. Such is the case of the 
NGC 4222 group, which may link the M group to Virgo at 
a mean distance ratio to Virgo of 1.2, and the NGC 4343 
group, around (12. 4h, 8°), similar to the B group of G99 
and introduced by de Vaucouleurs ( |1961| ) as its W group, 
and which may link the W group to the Virgo cluster at a 
mean distance ratio to Virgo of 1.5. A small group around 
NGC 4639 at (12. 7h, 14°) can also be identified with the 
eastern group of G99, although our group is less extended 
and less numerous; it could be linked to a remote group 
of galaxies around NGC 4746. 

In Table |, we propose a new nomenclature for all these 
groups, based on the name of its main galaxy, together 
with a list of the members we identified in this study. 
These lists are obviously not complete, because of the lack 
of distance measurements of some other well-known mem- 
bers. Table ^ lists other galaxies which have not been at- 
tributed to these groups nor to the main cluster, and are 



classified as background (class 4), infalling from behind 
(class 3), or from in front (class 1). 

After excluding all the galaxies which do not belong to 
the main cluster according to our classification, we are left 
with 71 spiral and 52 early-type galaxies. 67% of the origi- 
nal sample with distances are therefore considered as true 
members of the cluster; not surprisingly, the spiral sample 
exhibits a larger contamination (40%) than the early-type 
sample (20%). Among our 5 references, the proportion of 
true members varies from 56% (E00) and 61% (F98) to 
83% (T01). 

7. Comparison of Virgo distances 

Now that we have classified all the galaxies of our sam- 
ple into the different classes, we can return to the differ- 
ent distances adopted by the various authors to the Virgo 
cluster. 

We will first derive an average ratio of distances from 
a given reference to the predicted distance given by the 
Tolman-Bondi model and an adopted Virgo distance of 
15.4 Mpc. For the i th galaxy measured by reference j, we 
have: 



Tij ~ {d/d v )i x 15.4 1 ° J 

rj = (ni) (16) 

In the second test we restrict our sample to true Virgo 
galaxies (class 2) and compute the mean distance for 
each reference. The results will depend much less on our 
adopted Tolman-Bondi model but should confirm it. 

Table 4 gives the results. We can see that for G99 and 
T01, the agreement between the resulting mean Virgo dis- 
tances from both methods (Col. 3 and Col. 4, respectively) 
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Fig. 4. Histogram of the "Virgo derived distances" for spiral galaxies (left) and early-type galaxies (right) 



2000 - 




2000 



stooo - 




Fig. 5. Distribution of the recession velocities among the different classes for spiral galaxies (left) and early-type 
galaxies (right) 



and the value adopted by these authors (Col. 5) is satisfy- 
ing. On the contrary, for F98 the resulting mean Virgo dis- 
tances are consistent from both methods, but smaller than 
their adopted value: this is certainly due to inclusion of 
background galaxies into their sample; indeed, among the 
49 galaxies of their "fiducial sample" , 14 are classified by 
us as non-members of the Virgo cluster. For E00 the result 
is inverse and the mean error is large: it seems therefore 
that there is a discrepancy between thei r adopted Hub- 
ble constant (from Theureau et al. 
distance we derive from their data. 



1997]) and the Virgo 



For these two references, we adopt as the mean Virgo 
distance the average over the two determinations (Col. 3 
and Col. 4). This gives 18.9 and 23.2 Mpc for F98 and E00, 
respectively. If we apply to these distances the conversion 
factors adopted in Sect. 2.3.2 to reduce F98 and E00 to a 
mean Virgo distance of 16 Mpc, namely 0.832 and 0.678 
respectively, we now get values close to 16 Mpc (15.7 and 
15.8 Mpc, respectively). 

However, it is important to keep in mind that the clean 
samples of true members may st ill be affected by the in- 
completeness bias (Fouque et al. 1990 for the specific case 
of the Virgo cluster; Teerikorpi 1997 for a general discus- 
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Fig. 7. Distribution on the sky of the spiral galaxies (left) and early- type galaxies (right). Filled squares correspond 
to class 1, crosses to class 2, open squares to class 3 and open triangles to class 4 

Table 2. List of groups outside of the main Virgo cluster, with their mean distance ratio to Virgo and a list of 
members with known distance 



Name d/dv Location Members 



NGC 4168 


2.13 


M group 




2.07 


M halo 


NGC 4222 


1.24 


M infall 


NGC 4261 


2.11 


W group 




2.17 


W halo 


NGC 4343 


1.46 


W infall 


NGC 4639 


1.41 




NGC 4746 


1.92 





NGC 4168, NGC 4189, NGC 4193, NGC 4200, IC 769, IC 3061, IC 3099 

NGC 4067, NGC 4152, IC 3074 

NGC 4222, IC 3033, IC 3066, IC 3105, UGC 7249 

NGC 4180, NGC 4197, NGC 4215, NGC 4233, NGC 4235, NGC 4259, NGC 4260, 
NGC 4261, NGC 4273, NGC 4281, IC 3225, CGCG 42-1, CGCG 42-36 
IC 776, UGC 7579 

NGC 4252, NGC 4316, NGC 4318, NGC 4343, NGC 4353, NGC 4376, NGC 4390, 
NGC 4411A, NGC 4434, NGC 4451, NGC 4466, IC 3115, IC 3322A, UGC 7423 
NGC 4620, NGC 4633, NGC 4639, NGC 4651, IC 3742 
NGC 4746, UGC 8085, UGC 8114 



sion), which may lead to average cluster distances which 
are too short. To estimate the amount of this bias, we 
have calibrated once more the B-band Tully-Fishe r rela - 
tion using the 21 calibrators from Freedman et al. ( 2001 ), 
and the 51 true members of Virgo according to our classi- 
fication, which best match the calibrator properties (thus 
excluding peculiar, interacting, Hi-truncated galaxies, and 
restricting to morphological types between 2 and 8, incli- 
nations between 37 and 90 degrees, logVmax larger than 
1.7), with data extracted from the LEDA database. As 
usual, the slope of the Tully-Fisher relation is determined 
from the Virgo sample (more numerous) , and the intercept 
from the calibrators. The rms dispersion is similar for both 



samples (0.41 for calibrators, 0.53 for Virgo sample). The 
direct fit (appropriate to distance determinations) gives: 



M B = -6.28(±0.39)logV max - 6.31(±0.09) 



(17) 



The intercept for the Virgo sample, 24.96±0.11 gives a 
distance modulus of 31.27±0.14, or a distance of 18. Oil. 2 
Mpc, which is larger than our adopted distance from 
Cepheids. If we divide the sample of 51 Virgo galaxies 
equally into large and small logVmax, we get mean dis- 
tances of: 



H = 31.24 ±0.09, for log V max > 1.992 
H = 31.28 ± 0.12, for log V max < 1.992 



(18) 
(19) 
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Table 4. Mean ratio of galaxy distances to Tolman-Bondi predictions for dvirgo = 15.4 
Mpc, resulting mean Virgo distance, comparison with the mean distance restricted to 
true Virgo galaxies and to the distance adopted by each of our 5 references 



Reference (Galaxy type) 


Ratio 


dv 


dT c (N, a) 


jrcf 

d v 


G99 (E) 


1.017 ±0.031 


15.7 ±0.5 


15.3 ± 0.5 (36, 2.94) 


16.0 


T01 (E) 


1.084 ± 0.024 


16.7 ±0.4 


17.0 ± 0.4 (29, 1.99) 


17.0 ±0.3 


G99 (S) 


1.025 ±0.025 


15.8 ±0.4 


15.7 ±0.5 (41, 3.10) 


16.0 


F98 (S) 


1.221 ±0.027 


18.8 ± 0.4 


19.1 ± 0.5 (67, 4.27) 


21.5 ± 2.3 


E00 (S) 


1.487 ±0.053 


22.9 ± 0.8 


23.6 ± 1.1 (23, 5.51) 


21.8 



It can therefore safely be concluded that the incom- 
pleteness bias does not affect significantly our sample of 
true Virgo members (it would lead to a smaller mean dis- 
tance for galaxies with small rotation velocities, which are 
more biased). 

It can still be argued that the bias may be hidden, 
because the slope of the Tully-Fisher relation, when de- 
termined from the Virgo cluster sample, may be biased 
due to a possible incompleteness of the sample at small 
log Vmax. This is however not observed, as the direct slope 
determined from the calibrators is —5.24 ± 0.66, while 
an incompleteness bias would predict a steeper slope as 
compared to the slope determined from the Virgo sample 
(-6.28 ±0.39). 

The 0.3 mag difference between the Cepheid and Tully- 
Fisher distances is worrying. Tenants of the long distance 
scale will see evidence that the Tully-Fisher distance is the 
correct one, while most galaxies with Cepheid distances 
are infalling on the front side of the cluster. 



8. Conclusions 

In this paper, we have investigated how the relativistic 
Tolman-Bondi model as applied in E99 gives constraints 
on the Virgo cluster mass and distance, and allows one to 
disentangle its quite intricate structure. Distances to 183 
Virgo galaxies from 5 references have been used and aver- 
aged, together with HI deficiency parameters for spirals, 
to classify the galaxies into 4 different distance classes: an 
infalling component in front of Virgo, the Virgo cluster it- 
self, an infalling component behind the main cluster, and 
background groups. The main results of the present study 
are: 

— Among the 6 galaxies in the Virgo area with Cepheid 
distance measures, NGC 4639 is well-known to be dis- 
crepant. Our model nicely explains this apparent dis- 
crepancy, and the resulting mean distance to Virgo 
using all 6 calibrators is 15.4 ± 0.5 Mpc. However, 



the mean Tully-Fisher distance derived from 51 spi- 
ral galaxies classified as true members of the cluster 
in the present study, and calibrated using the same 
Cepheid distance system (Freedman et al. 2001 ) is 
larger (18.0 ± 1.2 Mpc). 

— The mass of the Virgo cluster derived from our model 
is large: M = 1.2 10 15 Mq, which corresponds to 1.7 
virial mass. Our adopted mass profile is steeper than 
the light distribution (anti-biasing). 

— Apart from the well-known background groups which 
once contaminated Virgo samples, we have identified a 
number of galaxies at intermediate distances between 
these groups and the cluster itself, which we interpret 
as filaments extracted from the background groups and 
falling into the cluster from behind. 

— We have not been able to securely identify galaxies in 
front of the cluster an d falli ng into it. This supports 
Sandage & Tammann ( 1976j ) who always argued that 
apparently foreground galaxies in the Virgo cluster di- 
rection were in fact cluster members. 

— Distances to the Virgo cluster adopted by F98 and 
E00 differ from what would be expected using the dis- 
tances predicted by our model and the mean ratio of 
published to predicted distances. They also differ from 
the mean of the published distances restricted to the 
true Virgo members. We also note that the distances 
published in E00 arc highly dispersed, even for true 
Virgo members. 
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